Features of pulsed synchronization of a 
systems with a tree-dimensional phase space. 

A. P. Kuznetsov 1 , N.V. Stankevich 2 and L.V. Turukina 1 

1 Institute of Radio Engineering and Electronics, Russian Academy of 

Sciences, Saratov Branch 
Zelenaya 38, Saratov, 410019, Russian Federation 
2 Saratov state University, Astrachanskaya, 83, Saratov, 410012, Russian 

Federation 

Abstract 

Features of synchronization picture in the system with the limit 
cycle embedded in a three-dimensional phase space are considered. 
By the example of Rossler system and Dmitriev - Kislov generator 
under the action of a periodic sequence of 5 - function it is shown, that 
synchronization picture significantly depends on the direction of pulse 
action. Features of synchronization tons appeared in these models are 
observed. 

1 Introduction 

Synchronization is a fundamental nonlinear phenomenon. It is presently the 
subject of intensive interest in many areas of a science pQ-jH]- The classical 
case of synchronization consists in an external periodic (usually harmonic) 
signal acting upon an autooscillating system with a stable limit cycle [3] , [6] , 
P - [12]. In this case, frequency locking and quasiperiodic regimes are possi- 
ble inside and outside the Arnold tongues, respectively, on the external signal 
frequency - amplitude plane. In phase space, these regimes are represented 
either by the stable torus or by the stable and saddle limit cycles appearing 
on this torus. However, there is a problems when external action has the 
form of short pulses. This external action can be represented as periodic 
sequence of 6— functions. Analysis of these problems is important both for 
the applications in many areas of science and for the theory of oscillations 
and nonlinear dynamics, where it is interesting to reveal possible specific 
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features of the synchronization picture. In series of the works [J3] - [19J this 
situations were studied within the framework of model system with a circle 
limit cycle under the action of the periodic sequence of 5 - function. For 
this model approximated one- dimensional circle map for the phase is con- 
structed. Authors was shown that synchronization picture both in a model 
system with a circle limit cycle and in a one-dimensional map had a number 
of features. This synchronization picture essentially different from the clas- 
sical as in sine circle map [6], [20J. These features have a certain degree of 
universality, since they refer to the investigation of a system in the vicinity 
of the Andronov - Hopf bifurcation. For example, in the papers [21], [22] 
similar results was obtained for the Van der Pol and Van der Pol - Duffing 
oscillators under the action of the periodic sequence of 5 - function. At the 
same time, in systems with two-dimensional phase space the external pulses 
can be directed along any direction in the plane, where the phase portrait 
of these systems represents a circle or the object close to circle. There are 
more complicated situation in the case then system with limit cycle embed- 
ded in tree-dimensional phase space [23], (21] • Let there is a plane in which 
limit cycle is the circle (which is usually the case near the Andronov - Hopf 
bifurcation). Then, external pulse acts in any direction in this plane lead 
to described above results. Different situation takes place when the external 
pulse acts in the direction perpendicular to this plane. Below, this situation 
is analyzed in application to the standard model systems in the nonlinear dy- 
namics. There are Rossler system, which is a artificial system, and Dmitriev 
- Kislov generator, which is described real electronic circuit. 

2 Pulsed synchronization of a Rossler system 

Let us first consider features of pulsed synchronization by the example of the 
Rossler system. We start with the case where the external action (a periodic 
sequences of ^-functions) acts in the plane in which the limit cycle of an 
autonomous Rossler system is situated. Then, the nonautonomous system 
can be described as 



where x, y, z are the dynamical variables; p, q, r are the system parameters; 
and A and T are the amplitude and period of the external action, respectively. 
Let us select p, q, r so as autonomous Rossler system has a single limit cycle of 
period 1, which is predominantly situated in the (x, y) plane and only a small 
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Figure 1: Attractor for the autonomous Rossler system. Attractor is plotted 
for p = 0.2, q = 0.1 and r = 1.5. 

part projects from this plane in the z axis direction (fig.l). Depend upon the 
parameters values altitude of the cycle "projection" in the z axis direction 
can be different. As we shall show further, the synchronization picture in 
nonautonomous system depends of these "projection" altitude too. 

Let us p = 0.2, q = 0.1 and r = 1.5. The altitude of the cycle "projection" 
in the z axis direction in this case is Az ~ 3.5. The period of oscillations 
in the autonomous system can be numerically evaluated as T = 5.4368309... 
Figure 2 shows the chart of dynamical regimes in system ([1]) and its fragments 
constructed using the results of computer simulation on the A — T parame- 
ters plane. In this and analogous charts to follow, the region painted white 
corresponds to the regime of period 1, bright-gray color corresponds to the 
regime of period 2 and above, black color refers to chaos and quasiperiodic 
regimes and dark-gray color indicate the region where the trajectory goes 
to infinity (the character of each regime is determined in the corresponding 
Poincare section). The chart in Fig. 2 exhibits a quite traditional pattern of 
the Arnold tongs. The increased fragments of this chart reveal the following 
features. For the large periods of the external action (in comparison to the 
period of autonomous system T = 5.43...), the pattern of synchronization 
is close to classical that corresponds to the simplest one-dimensional sine 
circle map (fig. 2c) [6J. The internal structure of synchronization tongs is 
traditional. One can see the period-doubling bifurcation and characteristic 
configurations " crossroad area " [20J. On the other hand, for external ac- 
tions with a period close to the eigenvalue, the main synchronization tong 




Figure 2: The chart of dynamical regimes and its fragments on the A — T 
parameters plane, constructed for nonautonomous Rossler system (pQ) with 
the parameters p = 0.2, q = 0.1 and r = 1.5. Region I corresponds to 
trajectories going to infinity. 

(albeit surrounded by a systems of tongs having the classical shape) exhibit 
a different internal structure (fig. 2b). This structure reveals a clear line of 
the secondary Andronov - Hopf bifurcation and the corresponding system of 
secondary synchronization tongs. 

Figure 3 shows the nonautonomous attractors of Rossler system (DO). 
Figures 3a and 3c shows various configurations of the period 3 attractors, 
and figures 3b and 3d shows the period 2 and 4 attractors accordingly. All 
figures concern to when the amplitude of external action is rather 

essential. If the amplitude is small, nonautonomous attractors will be dis- 
tinguished not strongly from an autonomous case. The external pulse not 
strongly changed a trajectory of autonomous system. Thus, the period of 
noautonomous system increased, but the attractor shape will change a little. 
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Figure 3: The portraits of nonautonomous attractors of Rossler system (TjQ). 
Attractors are plotted for the next pulsed amplitude and period: a) A = 4.4 
and T = 1.8; b) A = 1.8 and T = 10.9; c) A = 2 and T = 3.7; d) A = 2.5 
and T — 7.1. 

Now we change parameters of the system so, that the altitude of the 
cycle "projection" in the z axis direction increased. Let us p = 0.3, q = 0.1, 
r = 2.0. In this case Az w 6.7. The chart of dynamical regimes and its 
fragments on the A — T parameters plane, constructed for nonautonomous 
Rossler system (CQ) for this case are presented in Fig. 4. One can see, that in 
this well as in previous, synchronization tongs are observed. However, 

the region in which the trajectories going to infinity has increased. And the 
periodic and quasiperiodic regimes are realized at smaller values of amplitude 
of external action(fig. 4a). Also Arnold tongs were a little transformed (fig. 
4 b, c). They became narrower, their symmetry is broken, and they are a 
little bit inclined in the left side. 

Now let us consider how the synchronization picture in a nonautonomous 
Rossler system changes when the pulses act along the z axis, that is perpen- 
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Figure 4: The chart of dynamical regimes and its fragments on the A — T 
parameters plane, constructed for nonautonomous Ressler system (pQ) with 
the parameters p = 0.3, q = 0.1 and r = 2.0. 

dicularly to the plane to which the limit cycle predominantly belong: 
x = -y-z, 

V = x + py, (2) 
z = q + (x — r)z + AJ28(t — nT), 

Again we shall start with when autonomous Rossler system has a 

single limit cycle of period 1 and altitude of the cycle "projection" in the 
z axis direction is Az ~ 3.5, i.e. p = 0.2, q = 0.1 and r = 1.5. Figure 
5 shows the chart od dynamical regimes and its magnified fragments for 
this case. A comparison between Figs. 2 and 5 shows substantial changes 
in the pattern of synchronization. In particular, the map of system in 
Fig. 2 exhibits, in addition to the system of synchronization tongs and the 
region of quasi-periodic regimes (situated in the bottom part of the chart), 
rather large regions of the regimes of period 1 and 2, which correspond to 
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Figure 5: The chart of dynamical regimes and its fragments on the A — T 
parameters plane, constructed for nonautonomous Rossler system (J2D with 
the parameters p = 0.2, q = 0.1 and r = 1.5. 

higher values of the amplitude of external action. The lower boundary of the 
regimes of period 1 represents the line of the Neimark - Sucker bifurcation. 
The Arnold tongs touching this line have sharp vertices both at the base 
and the top (Fig. 5b). There is a small system of tongs forming open rings 
in the right-hand part of the chart. The other tongs on the right from the 
main one appear as bands. For large amplitudes these tongs transform into 
large regions corresponding to the regimes of various periods (Figs. 5a and 
5c). Note that in the case of small periods of external action the threshold 
from which the trajectories go to infinity shifts to zero. This behavior also 
differs from that observed in Fig. 1. The nonautonomous periods 1, 2, 3 and 
5 attractors of Rossler system (T5]) are presented in Figure 6 which should be 
compared to Fig. 3. 

Now let us changed parameters of the system so, that the altitude of the 
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Figure 6: The portraits of nonautonomous attractors of Rossler system (j2J). 
Attractors are plotted for the next pulsed amplitude and period: a) A = 5 
and T = 2.85; b) A = 2.4 and T = 2.85; c) A = 0.5 and T = 2.85; d) A = 1 
and T = 4.5. 

cycle "projection" in the z axis direction increased and Az rj 6.7. Figure 
7 shows the corresponding chart of dynamical regimes and its magnified 
fragments. 

A comparison between Figs. 4 and 7 shows that the region in which the 
trajectories go to infinity has increased and region of the of periodic and 
quasi-periodic regimes has decreased. This is especially well visible at large 
values of periods of external action. Besides synchronization tongs become 
narrower and are a inclined in the left side. They internal structure becomes 
complicated. The region of small values of period of external action (T < 6.0) 
is exception. Here the region of periodic and quasi-periodic regimes increased. 
And there are transition to chaos through period-doubling bifurcation inside 
synchronization tongs. 
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Figure 7: The chart of dynamical regimes and its fragments on the A — T 
parameters plane, constructed for nonautonomous Rossler system (|2j) with 
the parameters p = 0.3, q = 0.1 and r = 2.0. 

3 Pulsed synchronization of a Dmitriev - Kislov 
generator 

Now let us consider features of pulsed synchronization by the example of the 
Dmitriev - Kislov generator [20] . |24j : 

x + Dx = Mzexp(—z 2 ), 

y = x-z, (3) 

where x, y, z are the dynamical variables, D, M, Q are the system param- 
eters. Dynamics of this system is studied adequately. Thus we selected the 
system's ([3]) parameters so that in phase space there is the period-1 stable 
limit cycle. Let us fix D = 6, M = 5.5 and Q = 10. The tree-dimensional 
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attractor for system and its projections are presented in Fig. 8. One can 




Figure 8: Attractor and its projections for the Dmitriev - Kislov generator 
([3]). Attractor and projections are plotted for D = 6, M = 5.5 and Q = 10. 

see that projection of the attractor to the (y, z) plane is circle (Fig. 8 b). 
Two other projections is different form the circle, essentially, projection to 
the (x, y) plane. It is looks like the figure-of-eight (Fig. 8 a). 

Let us add into Dmitriev - Kislov generator external pulses in the form 
AJ2$(t — nT). Here A is the pulsed amplitude and T is the pulsed period. 
Note, that external pulses add each equations by turns. Thus, let us consider 
next situations. First case, when pulses act along the x axis, that is external 
action add to the fist equation. Second case, when pulses act along the y 
axis, that is external action add to the second equation. Third case, when 
pulses act along the z axis, that is external action add to the third equation. 
Figure 9 shows the charts of dynamical regimes on the pulsed amplitude - 
period plane. These charts are plotted for all cases referred above. Chart of 
dynamical regimes is the plane of parameters, at which regions of periodic 
behavior with different periods are shown with different colors. In this and 
analogous charts to follow, the region painted white corresponds to the regime 
of period 1, bright-gray color corresponds to the regime of period 2 and above, 
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Figure 9: Charts of dynamical regimes of the nonautonomous Dmitriev - 
Kislov generator on the A - T plane. Charts are plotted for the cases, when 
pulses act along the x axis (fig. a), the y axis (fig.b) and the z axis (fig.c). 
Systems parameters are fixed at D = 6, M = 5.5 and Q = 10. 



black color refers to chaos and quasiperiodic regimes. As can be seen from 
Figs. 2 and 3, the case when pulses act along the y axis is most interesting. In 
this case synchronization picture differ both from the classical typical for the 
sine circle map and from the synchronization picture typical for system with 
a circle limit cycle under the action of the periodic sequence of 5 - function. 
We can explains this fact as follows. In the case then external action add to 
the second equation of the system (to the dynamical variable y) pulses 
will be crossed both loops of the figure-of-eight (Fig. 8 a). In the other 
cases pulses are crossed either circle or one loop of the figure-of-eight. As 
a result, synchronization picture in the last two cases is the same like in 
the model with circle limit cycle under the action of the periodic sequence 
of 8 - functions (Figs. 3a and 3b). Thus, let us discuss synchronization in 



11 



nonautonomous Dmitriev - Kislov generator in the case then external pulses 
act along the y axis more detailed. 

The fragments of the chart of dynamical regimes of nonautonomous at- 
tractors for the Dmitriev - Kislov generator are presented in the Figs. 10 
and 11. In the charts one can see synchronization tongs corresponded to 




Figure 10: Fragment of the chart of dynamical regimes for the nonau- 
tonomous Dmitriev - Kislov generator ([3]) from the Fig.9b. The portraits 
of nonautonomous attractors are shown in this figure too. Attractors are 
plotted for the next pulsed amplitude and period: a) A = 3 and T = 0.8; b) 
A = 6 and T = 3.3; c) A = 3.7 and T = 3.13; d) A = 2 and T = 3.14; e) 
A = 1.5 and T = 4.1; f) A = 3 and T = 3. 

a periods multiple to the eigenvalue (the natural period of the system (jHJ) 
T ~ 6.22). Region of a quasiperiodic regimes are surrounded these tongs. 
Higher-order synchronization tongs are observed in this region. Synchroniza- 
tion tongs corresponded to subharmonic resonance are located in the region 
of lesser periods of external action (comparison to the natural period). These 
synchronization tong looks like strongly stretch up "drops". But note that, 
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Figure 11: Fragments of the chart of dynamical regimes for the nonau- 
tonomous Dmitriev - Kislov generator ([3]) from the Fig. 9 b. 

synchronization tong with rotation number | (period 2) differ from the oth- 
ers synchronization tongs in this region. This tong is touched to A = 
axis. First, it is converged in the middle part by forming "bridge". Then 
it is widen in the top part. Besides, in the top part of the chart this syn- 
chronization tong comes to the period 1 region (Fig. 4). Thus, on the one 
hand it is synchronization tong with rotation number |. On the other hand 
it is the period 2 region, which is originated from the regime of period 1 
being at large values of the pulsed amplitude as a result of period-doubling 
bifurcation. Nonautonomous attractors plotted for the Dmitriev - Kislov 
generator is verified this conclusion (Fig.4a-f). Thus, attractor plotted in 
the upper part of the period 2 synchronization region (Fig. 10 b) looks like 
period 2 attractor originated from the period 1 attractor (Fig. 10 4a) as 
a result of period-doubling bifurcation. But, attractor plotted in the lower 
part of the period 2 synchronization region (Fig. 10 c) looks different. It is 
corresponded to regime of synchronization with rotation number |. Transi- 
tion between these to types of the nonautonomous attractors observed under 
decrease of pulsed amplitude. This transition consist in the next. One at- 
tractor loop (see Fig. 10 b) is decreased throughout the height. While in the 
"bridge" it is situated in the same plane that action (Fig. 10 c). After that 
loop is lower oneself under this plane (Fig. 10 d). It is important note, that 
this process is accompanied by decrease of pulsed amplitude. These features 
are not observed for all others synchronization tongs in this region (Fig. 10). 

Synchronization tongs of different periods and quasiperiodic regimes are 
observed on the right of main tong (Fig. 11) However, these tongs are smaller 
by the amplitude and had more complicated structure. Inside these tongs 
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there are regions of different period and chaos, i.e. transition to chaos through 
period-doubling cascade. Moreover, there are other regions of the complex 
dynamics in the top part of the charts presented in Fig. 11. Thus, one 
can conclude that in the nonautonomous system ([3]) there are periodic and 
quasiperiodic regimes at the small values of the pulsed amplitude and periodic 
regimes and chaos at the large values. 

4 Conclusion. 

Thus by the example of Rossler system and Dmitriev - Kislov generator the 
problem of synchronization by the periodic sequence of 5 - function in the 
system with a three-dimensional phase space are considered. Different meth- 
ods of the adding of external action to the system are studied. It is showed 
that the synchronization in the system with limit cycle embedded in a tree- 
dimensional phase space significantly depends on the direction of pulse action 
and on the attractor's form. If external pulses act in the plane, where the 
projection of phase portrait of autonomous system represents a circle (or the 
object close to a circle), then synchronization picture in nonautonomous sys- 
tem will be closed to the synchronization picture typical for two-dimensional 
system with a circle limit cycle under the acting of the periodic sequence of 

5 - function. However, if projection of phase portrait essentially differ from 
the circle or if the external pulses act in the direction perpendicular to the 
plane of the limit cycle, then synchronization picture will be different. In the 
present paper it is showed that this synchronization picture will be differed 
both from the classical typical for the sinus-circle map and from the synchro- 
nization picture typical for two-dimensional system with a circle limit cycle 
under the action of the periodic sequence of 5 - function. 

This study was supported by the Russian Foundation for Basic Research 
(project no. 06-02-16773) and CRDF BRHE (grant REC-006 Y2-P-06-13.) 
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